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ABSTRACT 

In this letter, we introduce a new method of image stacking to directly study 
the undetected but possible 7-ray point sources. Applying the method to the 
Australia Telescope 20 GHz Survey (AT20G) sources which have not been de- 
tected by Large Area Telescope (LAT) on Fermi Gamma-ray Space Telescope 
(Fermi) , we find that the sources contribute (10.5±1.1)% and (4.3±0.9) % of 
the extragalactic gamma-ray background (EGB) and have a very soft spectrum 
with the photon indexes of 3.09±0.23 and 2.61±0.26, in the 1-3 and 3-300 GeV 
energy ranges. In the 0.1-1 GeV range, they probably contribute more large fac- 
tion to the EGB, but it is not quite sure. It maybe not appropriate to assume 
that the undetected sources have the similar property to the detected sources. 

Subject headings: gamma rays: diffuse background — methods: statistical — quasars: 
general 



Introduction 



The EGB was first detected by the SAS-2 mission (IFichtel et al.l 119751 ). and its spec- 



trum was rneasure d with good accuracy by Fermi (also called isotropic diffuse background. 



Abdo et al.ll2010dl ). It has been found to be consistent with a featureless power law with a 



^National Astronomical Observatories, Yunnan Observatory, Chinese Academy of Sciences, Kunming 
650011, China 

^Key Laboratory for the Structure and Evolution of Celestial Objects, Chinese Academy of Sciences, 
Kunming 650011, China 

•^Graduate School, Chinese Academy of Sciences, Beijing, P.R. China 



-2- 



photon index of ~2.4 in the 0.2-100 GeV energy range and an integrated flux (E>100MeV) 
of 1.03x10"^ ph ' 
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The origin of the EGB is one of the fundamental unsolved probl ems in astrophysics, 
and it has been a subject of study for a long time (see iKneiskd l2008l . for a review). The 
EGB could originate from either truly diffuse processes or from unresolved point sources. 
Truly d iffuse emission can arise frora nuraerous processes such as the annihilation of dark 
matter (lAhn er al.ll2007l : ICuoco et al.ll2010t iBelikov fc Hooperll2010l ) , particle acceleration by 
intergalactic shocks produced during large scale structure formation ( iGabici fc Blasii J 120031 ) 
etc. 

Blazars (including BL Lac objects, flat spectrum radio quasars, or unidentiffed flat 
spectrum radio sources) represent the most numerous population detected by the Ener- 
getic Gamma Ray Ex periment Te lescope (EGRET) on Compton Gamma Ray Observatory 
( iHartman et al.lll999l ) and Fermi ( lAbdo et al.ll2010al ). Therefore, the blazars which have not 
been detected by the EGRET or LAT are the most likely candidates for the origin of the 
bulk of the EGB emission. Many authors have studied the luminosity function of blazars 
and showed th at the contribution of blazars to the EGRET EGB could be in the range from 



20 % to 100 % fStecker fc Sa 



2008: Kneiske &: Mannheim 



amon 



2008 



1996 



Narumoto fc Totani 112006: Dermer 2007: Cao fc Bai 



Inoue &: Totanill2009[ ). Nevertheless, starburst galaxy and 



non-bl azar radio loud active galactic nuclei can also contribute a fairly large fraction of t he 
EGB (IThompson et al.l 120071 : iBhattacharya fc Sreekumarll2009l : iBhattacharya et al.ll2009[ ). 



Recently, lAbdo et al.l ( 1201 Obi ) built a source count distribution at GeV energy and 
yielded that point sources which had not been detected by the LAT can contribute 23 % 
of the EGB. At the fluxes currently reached by the LAT, they ruled out the hypothesis that 
point-like sources (i.e. blazars) produce a large fraction of the EGB. 

However, if the property of undetected sources is not similar to the detected sources, 
these conclusions maybe not correct. Therefore, we apply an image stacking method to di- 
rectly study the undetected point sources. For a sample of possible 7-ray point so urces which 



have not been detected by the Fermi due to their faint fluxes or soft sp ectra (lAbdo et al. 



2010d) . we can stack a large number of them to improve the statistics (jAndo &: Kusenko 



2010h . If their fluxes are not too faint, we can derive their mean flux and photon index by 
Maximum Likelihood (ML) method. 
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2. Sample 



The AT20gE| (iMurphy et al.ll2010[ ) is the largest catalog of high frequency radio sources 
and contains 5890 sources with the flux at 20 GHz exceeding 40 mJy in the whole sky south 
of declination 0°. In the south sky, about 60 % (230 sources) of Fermi 1-year LAT AGN 
Catalogue (ILAC) sources are associated with the AT20G sources ( iGhirlanda et al.ll2010al ). 
Through studying the co rrelation between the 7-ray and radio flux density of AT20G sources, 
Ghirlanda et al.l ( 1201 Obi ) yielded that AT20G sources not detected by the LAT can contribute 
17% of the EGB. Therefore, we apply our methods to study them firstly. 

We exclude the sources identified as Galactic Hn regions. Galactic Planetary Nebulas 
and parts of the Magellanic Clouds in AT20G, th en the majority of sources (5808) we obtain 
are quasi-stellar objects (see lMurphy et al.ll2010l . but the optical properties of these objects 
have not been published). In order to minimize the influence of strong sources, we only 
use the sources that are at least 2° away from the nearest First Fermi-LAT catalog (IFGL) 
source and locate at high Galactic latitudes, |6| > 15°. Finally, we obtain 2900 sources to 
analyze their contribution to the EGB. 



Method 



The photon^ we used in our analysis are same as that used by Abdo et al. (2010c) to 
construct the IFGL, but ours in the 1-300 GeV ener gy range have small 68% containmen t 
radius (better than 1°) and little source confusion (see lAtwood et al.ll2009l : lAbdo et al.ll2009f ). 
In this procedure, the tools of gtselect and gtmktim^ are used. 

For stacking the images of all sources, we collect all photons that are at most 1° away 
from any source of our sample and then record their energies {Ei, in units of GeV) and 
angular distances {9i, in units of deg) between the photon and the source. The overlaping 
of sources have little influence on our method because these sources are very faint and can 
be regard as parts of diffuse background source, especially in stacked image. 

After that, we apply a ML method to derive the flux and photon index of the stacked 
point source. For simplicity, in our model there are only two sources, e.g. the diffuse 



^It is available online through Vizier (http://vizier.u-strasbg.frl 
^http://fermi. gsfc.nasa.gov/cgi-bin/ssc/LAT/LATDataQuery.cgi 
These and other tools we used in next are accessible at 



http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.htinl 
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background source and stacked point source. We assume that they all have power-law 

spectra, and the photon indexes arc 71 and 72, respectively. The fluxes density are fi (in 
units of [ph cm~^ s~^ GeV~^ deg~^]) and /2 (in units of [ph cm~^ s~^ GeV~^]), respectively. 
The emission can be described by 

where dN{9, E) is the photon number in the ranges of {9-6+dd) and (E-E+dE), PSF is the 
point spread function (in the units of deg~^), exposure is the integral of effective area over 
time (in units of [cm^ s]). PSF and exposure are derived from the tool of gtpsf. The emission 
must meet the relationship of 

[ C^^'^T^^ ^' ^ /2(Y^)-"-P5F(^,i?)]exposm-e(i?)27r^d^di? = iVo, (2) 

where Nq is the total number of photons in the E1-E2 energy range. Therefore, there are 
only three free parameters. In the practical calculation, we use 71, 72 and M. M is the 
number of photons contributed by the stacked source. Then /i and /2 can be described by 
71, 72 and M. 

The probabihty for a photon with {9i,Ei) is 



Pi 



exposure{Ei)27r9i , Ei Ei 



t^^^lG^^'"' + /2(Y^)-'^^^5i^(^.,^0]- (3) 



The likelihood is the probability of the observed data for a specific model. For our case, it 
is defined as: 

No 

L = UPr (4) 

i=l 

The logarithm of the hkehhood is much conveniently calculated 

i=l i=l ^ 

Because the last term is model independent, it is not useful for ML method. Neglecting the 
last term, we get 

InL = J]ln[A(^)-^^ + f2{Y^r^PSF{e„E,)]. (6) 

i=l 

We maximize numerically L to obtain the most probable parameters (/i, /2 and 71, 72) of 
these sources. 
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We use the likelihood ratio to test the hypothesis. The point source "test statistic" (TS) 
is defined as 

T5 = -2(lnLo-lnLi), (7) 

where Lq and Li is the likelihood without and with point source. The detected significance 
of a point source is approximately y/TSa (see iMattox et al.lll996l ). 



4. Result & Discussion 

The stacked source is estimated to has a photon index of 2.81 and integrated fiux of 
1.07x10"^ ph cm s . The TS is 129, corresponding to a significance of ~11(t. The mean 
fiux of these sources is 3.69xl0^^^ph cm~^ s~^, it is fainter than the faintest IFGL source 
by a factor of 10. It can contribute 8.4% of the EGB in the 1-300 GeV energy range. We 
also apply our method to a subsample of fiat spectrum radio sources (i.e. a(5-20GHz) <0.5, 
with Fi, oc z/~", 1780 sources). Its photon index is 2.79, only slightly harder than the former. 
Its mean fiux is 3.79x10"^^ ph cm~^ s~^, and the TS is 88. This subsample has not distinct 
characteristic from the other sources in 7-ray energy range. 

In order to test a more complicated spectral shape of the stacked source, we analyze the 
spectrum in the 1-10 GeV energy range . We expect that the spectrum would be harder in 
this energy range. However, the estimated photon index is 3.01. Therefore, we analyze the 
spectrum in 2-10 GeV, 3-10 GeV energy range, respectively. The results are summarized in 
Table 1. We find that the spectrum is very soft in 1-3 GeV energy range and becomes harder 
above 3 GeV. It is indicated that two types of sources exist, in which one with softer and 
another with harder spectrum in GeV range. The former will dominate in lower energies, 
and latter in higher energies. Therefore, the spectrum of stacked source shows very soft 
in the 1-3 GeV energy range. We will study this further if the optical properties of these 
objects can be obtained. 

Finally we obtain the properties of the spectrum as follows. In the 1-3 GeV, the photon 
index is 3.12, the mean integrated fiux is 3.89 xlO^^^ ph cm~^ s~^, and the TS is 92. Obvi- 
ously the fiux is larger than that in the 1-300 GeV energy range. It could be caused by that 
the spectrum in the 1-300 GeV is not well fitted with a single power-law. In the 3-300 GeV, 
the photon index is 2.66, the mean integrated fiux is 3.72 xlO~^^ph cm~^ s~^, and the TS 
is 38. 



An decrease in In L of 0.5 from its maximum value corresponds to the 68% confidence 
(1 a) region for each parameter (see iMattox et al.lll996l ). We use this variance to estimate 
the error of each parameter. In three parameters (71, 72 and M), we take two ones to be 
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the values with maximal likelihood, and allow third one to change around its best value, we 
then test the deviation of InL from its maximum value shown in figure 1. The 1 a errors of 
72 are 0.25 and 0.22, the ones of M are 270 and 53, in 1-3 and 3-300 GeV respectively. The 
1 a relative errors of fluxes are 10.5 % and 17.3 %, in 1-3 and 3-300 GeV. 

In order to verify the effectiveness and accuracy of our method, we do the Monte Carlo 
(MC) simulations using the tool gtohssim. The simulating time is 26 Ms, equaling to the 
time of real data we used. We simulate the Galactic and isotropic diffuse backgrounds using 
the models (e. g. gllJem_v02.fit, isotropicJem_v02.txt) recommended by the LAT team, in 
which 3913 sources are generated each time, but only 2900 sources isotropically distribute 
on the sky with |6| > 15°. We complete one thousand MC simulations in 1-3 and 3-300 GeV 
energy range using the obtained parameters. The diffuse source is simulated only once due 
to long run time, but its effect on the result is not remarkable because the source is random 
distribution and its photons are various. 

The distributions of the photon index, flux and TS for different energy ranges are shown 
in figure 2. They are compatible with Gaussian distributions. In the 1-3 and 3-300 GeV, 
the mean fluxes are 4.30 and 0.364 (in the unit of [10^^^ ph cm~^ their relative errors 

are 10.3% and 20.7%; the photon indexes are 3.15 and 2.71 with the errors of 0.23 and 0.26. 
The errors estimated here are similar to that found in the fourth paragraph. Comparing 
the input parameters, we find that the systematic errors occur, especially for the flux in the 
1-3 GeV energy range. They could be caused by that diffuse background source can not be 
described by a single power-law spectrum. 

Because the MC method can obtain the systematic errors, we use this method to correct 
our results as follows: in the 1-3 and 3-300 GeV, the fluxes are 3.48±0.36 and 0.380±0.080 
(in unit of [10~^^ph cm~^ s~^]), and the photon indexes are 3.09±0.23 and 2.61± 0.26 
respectively, while the contribution to th e EGB is (10.5ib l .l) % a nd (4.3±0.9) %, which are 



much smaller than the result (17%) of iGhirlanda et al.l ( l2010bl ). If the soft spectrum in 



1-3 GeV is caused by the spectral broken of some sources, the photon index would not be 
extrapolated to lower energy range. However, as long as the spectrum of stacked source is 
not harder than the EGB, the contribution to the EGB in 0.1-1 GeV will be larger t han that 



in 1-3 GeV. Our result is compatible with the result (23%) of lAbdo et al.l ( l2010bl ) because 
other point sources could contribute to the EGB. 

In this letter, we introduce a new method of images stacking to directly study the 
contribution of undetected point sources to the EGB. Our method is more direct than the 
methods used by many authors. Those methods involve the 7-ray luminosity of undetected 
sources which is estimated through the properties of a few detected sources. They include 
many uncertainties and lead the result to be questionable validity. 
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Applying our method, we find that the undetected sources in AT20G can contribute 
(10.5±1.1)% and (4.3±0.9) % to the EGB in the 1-3 and 3-300 GeV energy range respec- 
tively. Their 7-ray spectrum is very soft, implying that the emissive property is different for 
undetected and detected sources. 

Applying our method to estimate the contribution of all point sources to the EGB, we 
need a complete sample of possible 7-ray point sources which is not easily constructed. In 
this letter, we only estimate the contribution of AG20G to the EGB. We will study more 
samples of possible 7-ray point sources in the future. 
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Table 1: The results in different energy ranges. 

Energy Ranges Photon Indexes Mean Fluxes TS 



(GeV) (10-^^ ph cm-^ s-^) 

1-300 2^81 30 129 

1- 10 3.02 38.4 108 

2- 10 2.95 7.49 37 

3- 10 2.64 2.84 19 
1-3 3.12 38.7 92 

3-300 2.66 3.72 38 




Fig. 1. — The change of hkehhood from its maximum with 72 (top) or M (bottom) around 
their best value when other parameters are fixed on their highest hkehhood value. Right 
and left panels are for 1-3 GeV and 3-300 GeV bands, respectively. 
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Fig. 2. — Distributions of photon indexes (top), mean fluxes (middle) and TS (bottom). 
Right and left panels are for 1-3 GeV and 3-300 GeV bands, respectively. The distribu- 
tion can be represented by Gaussian functions (dashed line) with central values /i=3.15 
(2.71), 4.30 (0.364) xlO"", 110 (37), respectively, standard deviations a=0.23 (0.26), 4.44 
(0.755) xlO^i^ 21 (13), in the 1-3 (3-300) GeV. The input photon indexes are 3.12, 2.66, 
fluxes are 3.87, 0.372 (in unit of [lO^^^ph cm~^ s~^]), respectively, which are indicated by 
vertical line 



